arXiv:1501.07004v2 [cond-mat.quant-gas] 8 Oct 2015 


Spin-orbit coupled repulsive Fermi atoms in a 
one-dimensional optical lattice 

Xiaofan Zhou, Kuang Zhang, Junjun Liang, Gang Chen and 
Suotang Jia 

State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of 
Laser Spectroscopy, Shanxi University, Taiyuan 030006, China 

E-mail: liangjj@sxu.edu.cn, chengang971@163.com 
9 October 2015 

Abstract. Motivated by recent experimental development, we investigate spin-orbit 
coupled repulsive Fermi atoms in a one-dimensional optical lattice. Using the density- 
matrix renormalization group method, we calculate momentum distribution function, 
gap, and spin-correlation function to reveal rich ground-state properties. We find that 
spin-orbit coupling (SOC) can generate unconventional momentum distribution, which 
depends crucially on the filling. We call the corresponding phase with zero gap the 
SOC-induced metallic phase. We also show that SOC can drive the system from the 
antiferromagnetic to ferromagnetic Mott insulators with spin rotating. As a result, a 
second-order quantum phase transition between the spin-rotating ferromagnetic Mott 
insulator and the SOC-induced metallic phase is predicted at the strong SOC. Here 
the spin rotating means that the spin orientations of the nearest-neighbor sites are not 
parallel or antiparallel, i.e., they have an intersection angle 9 £ (0, n). Finally, we show 
that the momentum /c pea k, at which peak of the spin-structure factor appears, can also 
be affected dramatically by SOC. The analytical expression of this momentum with 
respect to the SOC strength is also derived. It suggests that the predicted spin-rotating 
ferromagnetic (fc pea k < 7r/2) and antiferromagnetic (7r/2 < fc pea k < 7r) correlations can 
be detected experimentally by measuring the SOC-dependent spin-structure factor via 
the time-of-flight imaging. 
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1. Introduction 

Ultracold Fermi atoms in optical lattices have attracted considerable interest both 
experimentally [UISlISlSliniEllTllHllQl [lOl HU H21 HSJ [1^1 [T5], |T6] and theoretically 
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because these setups are powerful platforms to simulate rich physics of strongly- 
correlated materials [42] One of advantages of this system is that the spatial 

geometry of optical lattices can be well controlled. Especially, using a strong 
harmonic transverse confinement, one-dimensional (ID) optical lattices have been 
achieved experimentally W On the other hand, the relative parameters have high 
controllability, and moreover, can reach the regimes that cannot be accessible in the 
conventional condensed-matter physics. For example, the two-body interaction between 
Fermi atoms can be tuned by a magnetic-field-dependent Feshbach resonant technique 
i ni . and thus ranges from the positive (repulsive) to the negative (attractive). For 
the on-site repulsive interaction, a well-known second-order quantum phase transition 
between an antiferromagnetic Mott insulator and a metallic phase can emerge [45] . 

Another important breakthrough in recent experiments of ultracold Fermi atoms 
is to successfully create a synthetic spin-orbit coupling (SOC), with equal Rashba and 
Dresselhaus strengths, by a pair of counter-propagating Raman lasers ma im m S9]- 
Indeed, SOC describes interaction between the spin and orbit degrees of freedom of a 
particle. In contrast to the typical property of solid state materials that the intrinsic 
SOC strength is generally smaller than the Fermi velocity of electrons, this synthetic 
SOC strength realized can reach the same order as (or even larger than) the Fermi 
velocity of atoms, and moreover, can also be tuned in a wide range [50, 5T]. Recent 
theory has revealed that SOC can generate exotic superfluids, including topological 
Bardeen-Cooper-Schrieffer (52, 53] [51, 55] 150, 53 EH, 59] and topological Fulde-Ferrcll- 
Larkin-Ovchinnikov [60l f6TI 62j phases, for the attractive Fermi atoms. The fundamental 
picture for achieving these nontrivial topological superfluids is that SOC, Zeeman held, 
and s-wave interaction can induce triplet p-wave pairing [63i [641 [65] . In parallel with 
the attractive case, it is natural to ask what novel physics can occur in the repulsive 
Fermi atoms driven by the synthetic SOC [66] E3 1 68] . 

Inspired by the above experimental developments and theoretical considerations, 
here we investigate spin-orbit coupled repulsive Fermi atoms in a ID optical lattice. 
Recently, spin-orbit coupled Bose-Einstcin condensates in the ID optical lattice has 
been prepared experimentally [69]. Using a similar technique, the system considered 
could also be achieved in the near future. Physically, the spin-orbit coupled Fermi 
atoms in the optical lattices have two characteristics. One is that SOC can make Fermi 
atoms hop between the nearest-neighbor sites with spin flipping; see the Hamiltonian 
(|6i) in the following. Moreover, it has a strong competition with the on-site repulsive 
interaction, and especially, with the conventional spin-independent hopping; see the 
Hamiltonian ([4| in the following. The other is that in the presence of SOC, for different 
chemical potentials, the fillings are quite different; see Fig. [0 

Notice that in ID quantum fluctuation becomes significant, and the mean-held 
results are, in principle, unreliable m Here we capture the required ground-state 
properties, including momentum distribution and spin-correlation, by using the density- 
matrix renormalization group (DMRG) method [711 [72] . which is a powerful numerical 
method to study lower-dimensional strong-correlated systems [73]. The main results 
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Figure 1 . Spin-orbit coupled energy bands in an optical lattice and the corresponding 
fillings for different chemical potentials. The blue arrows represent spin polarizations in 
different subbands. The red dot denotes the contact point of two subbands. fc m ; n and 
fcmax are the minimum and maximum of subbands, and are determined by Eqs. m 
and m, respectively, if the on-site repulsive interaction is not taken into account, 
(a) and (b) show the fillings for the small chemical potentials. In these cases, only 
the lower subband (the black shadow) is partly occupied. In (a), there are four Fermi 
points and no occupation around k = 0 occurs, while in (b), there are three Fermi 
points and the occupation at k = 0 emerges, (c) and (d) show the fillings for the large 
chemical potentials, in which the occupations in both the lower and upper (the red 
shadow) subbands emerge. In (c), the lower subband is still partly occupied, while in 
(d), the lower subband is fully occupied for a larger chemical potential. 


are given as follows. Section 2 is devoted to introducing our proposal and deriving a 
ID Fermi-Hubbard model with the synthetic SOC. Section 3 is devoted to addressing 
the generalized results without the on-site repulsive interaction. In this section, the 
unconventional momentum distribution, which depends crucially on the filling, are 
found. We call the corresponding phase the SOC-induced metallic phase. Section 4 
is devoted to discussing the results in the presence of the on-site repulsive interaction. 
By means of the spin-correlation function, we find that SOC can drive the system from 
a spin-rotating antiferromagnetic Mott insulator to a spin-rotating ferromagnetic Mott 
insulator. As a result, a second-order quantum phase transition between the spin- 
rotating ferromagnetic Mott insulator and the SOC-induced metallic phase is predicted 
at the strong SOC. Here the spin rotating means that the spin orientations of the 
nearest-neighbor sites are not parallel or antiparallel, i.e., they have an intersection 
angle 9 e (0,7r) (see Fig. |H] in the following). In the spin-rotating antiferromagnetic 
Mott insulator, tt/2 < 8 < tt, and the quasi-long-range spin correlation decays as a 
power law and changes the sign with a period 2 < T < 4, whereas for the spin-rotating 
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Figure 2. (a) Proposed experimental setup for realizing spin-orbit coupled repulsive 
Fermi atoms (FA) in a ID optical lattice, (b) The interaction between FA and a pair of 
counter-propagating Raman lasers, labeled respectively by Ri and R 2 . These Raman 
lasers create a momentum-sensitive coupling between two internal atomic states, i.e., 
a synthetic SOC with equal Rashba and Dresselhaus strengths, (c) Energy levels and 
their transitions induced by Raman lasers. 


ferromagnetic Mott insulator, 6 < tt/2, and the quasi-long-range spin correlation 
also decays as a power law but with the period T > 4. Finally, we show that the 
momentum fc pea k, at which peak of the spin-structure factor appears, can also be affected 
dramatically by SOC. The analytical expression of this momentum with respect to the 
SOC strength is also derived. It suggests that the predicted spin-rotating ferromagnetic 
(fcpeak < tt/2 ) and antiferromagnetic ( tt /2 < k pea k < 7 r) correlations can be detected 
experimentally by measuring the SOC-dependent spin-structure factor via the timc-of- 
flight imaging 0 . The discussions and conclusions are given in section 5. 

2. Model and Hamiltonian 

2.1. Proposed experimental setup 

Figure [2] shows our proposal that repulsive Fermi atoms in a ID optical lattice are 
driven by a pair of counter-propagating Raman lasers. For the specific experiments, 
3D optical lattice is first prepared by the interference of three pairs of counter- 
propagating laser beams [32103]. The corresponding periodic potential can be written 
as V3D = Vo cos 2 (k w x) + Vo cos 2 (k w y) + Vo cos 2 (k w z), where Vo is the lattice depth, 
k w = 2 tt/X w is the wave vector, and A™ is wavelength. By further using a strong 
harmonic transverse confinement V = m 0 ojj_r 2 / 2 in the 3D optical lattice, i.e., the 2D 
harmonic potential frequency co_ j_ is far larger than the trapping frequency oo z along 
the weakly-confining axis, the required ID optical lattice can be generated M ; see 
Fig. [2])a) . In such case, the two-body interaction between Fermi atoms is described 
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effectively by m 
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where C ~ 1.46, a± = (2/t/mocax) 1 / 2 , 030 is the 3D s-wave scattering length, and m 0 
is the atomic mass. I 11 addition, a pair of counter-propagating Raman lasers shown 
in Fig. [2jb) are used to create the required ID synthetic SOC, with equal Rashba and 
Dressclhaus strengths In this method, the corresponding two spin states are chosen 
as |t) = |9/2, —9/2) and |t) = |9/2, —7/2) for 40 K system {46], or |t) = |3/2, —3/2) and 
|t) = |3/2, —1/2) for 6 Li system [49]; see Fig. ]2])c). 


2.2. Hamiltonian 

The total dynamics illustrated by Fig.[2]is governed by the following ID Fermi-Hubbard 
model with the synthetic SOC [TJ, 72] : 


H = H t + H u + H soc 

(3) 

H t = -t (4x0+ la + H.C.), 

(4) 

1 , 0 = T,t 


H U = U Y n; t n 4 , 

1 

(5) 

Hsoc = a B4tO+4 - 4+O+it + H.C.), 

(6) 


1 


where c\ a and Cj a are the creation and annihilation operators, with spin a =t, 4, at lattice 
site l, riifj = c] a ci a is the number operator, t is the spin-independent hopping magnitude, 
U is the on-site repulsive interaction strength, and A is the SOC strength. Based on 
above proposal experimental setups, the relative parameters can be tuned independently. 
For example, the hopping magnitude t can be controlled by the intensities of lasers 
[42, [43], the ID on-site repulsive interaction strength U can be tuned by Feshbach 
resonance ra. and the SOC strength A can be driven through a fast and coherent 
modulation of Raman lasers [50, El]- As a consequence, for a proper optical lattice, the 
SOC strength A has the same order of the hopping magnitude t m 
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2.3. Momentum distribution and spin correlation 

In terms of the SOC-induced properties (see Introduction), here we mainly focus on 
momentum distribution and spin correlation, which can be measured experimentally 
by the timc-of-flight imaging [4, T5| 06, 74] . The momentum distribution functions for 
spin-up and spin-down atoms are written respectively as m 

«#) = j E ( c W) • (?) 

1,3 

= T E eit ^ n (4 C «) • (8) 

1,3 

For the Hamiltonian ((2]), the spin-up and spin-down atoms are equal. It means that 
rq-(fc) is the same as n±(k), and thus we only consider n^(k) in the following discussions. 
The spin-correlation function is dehned as [Tamms] 

= J E AMD , (9) 

l 

where r is a distance between different sites and sf = cLq| — cj^cq. The corresponding 
spin-structure factor is given by oimim] 

sw^E^W- (i°) 

1,3 

Since the spin-structure factor has the sum extending over all lattice sites l and j, it 
reflects spin correlation globally, and is thus used experimentally to detect the magnetic 
order [74] . 

In addition, we will perform the DMRG calculations, with open boundary condition, 
to calculate Eqs. (I7|)- (jl0|) . The basic energy scale is chosen as t — 1. In the detailed 
calculations, we retain 150 truncated states (which is sufficient) per DMRG block and 
20 sweeps with the maximum truncation error ~ 10 -5 [501 [81] . 

3. Without on-site repulsive interaction 

In order to better understand the fundamental physics induced by SOC, we first 
consider a simple case without the on-site repulsive interaction (U/t = 0), in which 
the Hamiltonian (|3]) reduces to //, = //, + // soc , he., 

Hi = —t ^2 (4aC l+ ia + H.c.) + A £(4 t D + 4 - cj^Q+it + H.c.). (11) 

Z,(7=t,t l 

In experiments, the ID noninteracting Fermi atoms can be realized by tuning the 3D 
s-wave scattering length a 3 o to its zero crossing. In such case, the ID s-wave scattering 
length a id —» oo, and the ID effective two-body interaction U(z) then becomes zero 
P [7, 43]; see Eqs. (JTj) and ([2]) . 
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Figure 3. (a)-(c) The momentum distribution functions n-f-(fc) for the different SOC 
strengths X/t, when the filling factors are chosen as (a) n = 0.3, (b) n = 1.0, and 
(c) n = 1.7. (d) The momentum distribution functions n^(k) for the different filling 
factors n, when the SOC strength is chosen as X/t = 1. In all subfigures, the on-site 
repulsive interaction strength and the lattice length are given by U = 0 and L = 100, 
respectively. 


3.1. Momentum distribution function 


In the absence of SOC (X/t = 0), the system have three known features. Firstly, the 
energy bands of the Hamiltonian H t are degenerate, and the system is located at the 
metallic phase. Secondly, all Fermi atoms occupy the degenerate energy bands from 
k miT1 = 0, where k min is the minimum of the energy bands. Moreover, there are two 
degenerate Fermi surfaces at the Fermi momentum kp = u7t/ 2 [82], with filling factor 
n = N/L , where N is the total atomic number. At last, the chemical potential cuts the 
degenerate energy bands, with the Fermi momenta ±kp. As a result, the momentum 
distribution function n^(k) has a plateau of 1, with sharp edges at k = ±kp] see the 
black solid curve of Fig. [3]]a). In the large-A: limit, this plateau of 1 disappears and two 
plateaus of 0 emerge, as expected. For the different chemical potentials, the fillings are 
similar, and the momentum distributions are thus similar; see the black solid curves of 
Figs. |3]]b) and[3]]c). 

In the presence of SOC (X/t ^ 0), the results are very interesting. In this case, the 
energy bands governed by the Hamiltonian Hi split into two nondegenerate subbands, 
whose minima and maxima are given respectively by 


k m in = ±2 arctan( 


-t + Vt 2 + A 2 , 


A 


( 12 ) 
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fcmax = ±2 arctan( 


t + VW^xs, 

X 


( 13 ) 


It is easy to find that fc mi n + k max = vr. For A = 0, k m \ n = 0, as expected. In the presence 
of SOC (A 7 ^ 0), there are generally four Fermi points, when the chemical potential 
cuts two nondegenerate subbands. If the chemical potential is identical to the critical 
chemical potential p that just cuts the contact point [see the blue curve in Fig. [l](b)], 
three Fermi points can emerge. More importantly, for the different chemical potentials, 
the fillings are quite different; see Fig. [0 These different Filings affect dramatically on 
momentum distributions. As examples, we plot, in Fig. [31 the momentum distribution 
functions n^{k) for the different SOC strengths X/t, when the filling factors are chosen 
as (a) n = 0.3, (b) n = 1, and (c) n = 1.7. 

It can be seen from Fig. |3](a) that for a smaller filling, any SOC leads to a new 
momentum distribution, in which the corresponding function n^(k) has two plateaus 
of 1/2, apart from the conventional plateau of 1. The physical explanation is given as 
follows. When the SOC strength is not strong enough (see, for example, X/t = 0.2), 
p > p, and thus both the lower and upper subbands are partly occupied; see Fig. [0(c). 
The occupation in the upper subband determines the plateau of 1, while the occupation 
in the lower subband governs two plateaus of 1/2. With the increasing of the SOC 
strength X/t, i.e., the chemical potential p decreases [83, El], the occupation in the 
upper subband becomes less, and thus the wide of the plateau of 1 becomes narrower. 
In particular, when X/t = 0.52, the plateau of 1 disappears, since in this case p — jl, 
and thus no occupation in the upper subband can be found. If further increasing the 
SOC strength X/t, i.e., p < p, a plateau of 0 emerges around k = 0 (not in the 
large- k limit). This is because there is no occupation around k = 0 when p < p; see 
Fig. d a). For the half filling (n = 1), we find that the momentum distribution function 
7Z|(/c) usually has two plateaus of 1 and 1/2 [see Fig. [31(b)], because p > p and both 
the lower and upper subbands are partly occupied. For a larger filling factor (see, for 
example, n = 1.7), the momentum distribution function rp(k) has two plateaus of 1 
and 1/2 when X/t = 2.0; see Fig. |3](c). In this case, p p and the lower subband is 
totally occupied with a part occupation in the upper subband; see Fig. d(d). This case 
cannot occur in the smaller filling factors. We call the corresponding phase, with above 
unconventional momentum distributions, the SOC-induced metallic phase. Finally, in 
order to see clearly the above evolution of momentum distribution with respect to the 
filling factor, we plot the momentum distribution functions rp/k) for the different filling 
factors n in Fig. [3](d). 

It should be noticed that in the case with a Zeeman field, but without SOC, the 
chemical potential cutting the Zeeman-split bands also results in more than two Fermi 
points. But the momentum distribution function n^(k) is quite different from that 
induced by SOC. In such case, the momentum distribution function n^(k) only has the 
plateaus of 1. Moreover, for the different chemical potentials, the fillings are similar, 
and thus the momentum distribution functions n^(k) are also similar. However, in our 
considered case with SOC, but without the Zeeman field, the momentum distribution 
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Figure 4. The spin correlations between (a) the nearest-neighbor sites, i. e -» ( s ? s ?+ 1). 
and (b) the next-nearest-neighbor sites, i.e., (sfsf +2 ), for the different SOC strengths 
X/t. (c) The spin-structure factors S(k) for the different SOC strengths X/t. In all 
subfigures, the filling factor, the on-site repulsive interaction strength, and the lattice 
length are given by n = 1, U = 0, and L = 100, respectively. 


function n^(k) has the new plateaus of 1 / 2 , apart from 1 , and moreover, the filling- 
dependent unconventional momentum distributions emerge. 

3.2. Spin-correlation function and spin-structure factor 

Since SOC can make Fermi atoms hop between the nearest-neighbor sites, with 
spin flipping [see Fig. Ha) and the Hamiltonian (| 6 j)], it has a competition with the 
conventional spin-independent hopping. This competition has a strong effect on spin 
distributions of different sites. To see this clearly, we consider the spin-correlation 
function in Eq. (jUJ). 

In Figs. Ha) andHT), we plot the spin correlations between the nearest-neighbor 
sites, i.e., (sfsf +1 ), and the next-nearest-neighbor sites, i.e., (sfsf +2 ), respectively. These 
two subfigures show clearly that in the absence of SOC (X/t = 0), (sfsf +1 ) < 0 and 
( s i s i+ 2 ) = 0; which means that this short-range spin correlation is negative. In the 
presence of SOC ( X/t 7 ^ 0), the spin correlation is still short range. Interestingly, with 
the increasing of the SOC strength X/t, the spin correlation varies from the negative to 
the positive at the critical point A c /t = 1. The physical reason will be illustrated in the 
following section. 

I 11 Fig. H c )> we pl°t tl ie spin-structure factors S(k) for the different SOC strengths 
X/t at the half filling (n — 1). In the absence of SOC (X/t = 0), the spin-structure factor 
S(k) evolves as a straight line from 0 to 0.25, and has a cusp at k = ± 7 r [73]. For a finite 
SOC strength X/t = 0.5, the spin-structure factor is 5(0) > 0 and 0 < S(n) < 0.25. 
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Figure 5. The momentum distribution functions rp(k) for the different on-site 
repulsive interaction strengths UIt (Left panel), the gap A as a function of U/t (Center 
panel), and the gaps A as functions of 1/L for the different U/t (Right panel) at the 
half fillings (n = 1). In all subfigures, the SOC strength is chosen as (a-c) X/t = 0 and 
(d-f) X/t = 1.0. In (a), (b), (d), and (e), the lattice length is given by L = 100. The 
scaling behaviors in (c) and (f) show that in the thermodynamical limit (L —> oo), 
A = 0 for U/t = 0. 


When X/t = 1, the spin-structure factor S(k) becomes a constant 0.125, since in this 
case no short-range spin correlation can be found. When X/t = 2, the spin-structure 
factor S(k) is similar to the case of X/t = 0.5. 

4. With on-site repulsive interaction 

We now consider the case with the on-site repulsive interaction ( U/t > 0). In the absence 
of SOC (X/t = 0), it has been demonstrated exactly that for the ID homogeneous 
Fermi-Hubbard model at the half filling (n — 1), the metallic phase occurs at U — 0. 
For U > 0, the system is always located at the antiferromagnetic Mott insulator 
in which the spin orientations of the nearest-neighbor sites are antiparallel [78] • It 
means that a second-order quantum phase transition between the metallic phase and 
the antiferromagnetic Mott insulator emerges at U — 0. In the following, we will 
show that SOC can drive the system from the antiferromagnetic to ferromagnetic Mott 
insulators with spin rotating (i.e., the spin orientations of the nearest-neighbor sites 
are not parallel or antiparallel), and predict a second-order quantum phase transition 
between the spin-rotating ferromagnetic Mott insulator and the SOC-induced metallic 
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Figure 6. (a) The momentum distribution functions n-j-(fc) for the different filling 
factors n. (b) The gap A as a function of the filling factor n. (c) The gap A as 
functions of 1/L for the different filling factors n. In all subfigures, the SOC strength 
and the on-site repulsive interaction strength are given by \/t = 1.0 and U/t = 10, 
respectively. In both (a) and (b), the lattice length is given by L = 100. The scaling 
behavior in (c) shows that in the thermodynamical limit (L —> oo), A = 0, except for 
the half filling (n = 1). 


phase at the half filling {n = 1 ). 

4-1. Momentum distribution function and gap 

Figures [5](a) and[5](d) show the momentum distribution functions n^{k) for the different 
SOC strengths X/t at the half filling (n = 1 ). In the absence of SOC (A /t = 0) and the 
on-site repulsive interaction strength ( U/t = 0 ), the momentum distribution function 
n^(k) has two sharp edges at k = ±/c F ; see the black solid curve in Fig. [5](a) . When 
increasing the on-site repulsive interaction strength U/t , the momentum distribution 
functions n^(k) become smoother and all sharp edges disappear. In order to see the 
relevant physics more clearly, we introduce the gap [35] 

A = tL + -p~, (14) 

where pU = E g {N + 1) — E g (N ) and p~ = E g (N) — E g (N — 1), with the ground-state 
energy E g . This gap reflects the difference between the energy required to add (/i + ) 
and remove (/i~) a Fermi atom from the ground state. In the thermodynamic limit 
(L —> oo), it has been demonstrated rigorously that A = 0 for U/t = 0 and A 7 ^ 0 for 
U/t > 0. Moreover, A = 0 corresponds to the metallic phase and A / 0 corresponds to 
the Mott insulator [45]. In our numerical results, due to f ini te-size effects, the gap is not 
absolute zero when U/t = 0; see Fig. [5](b). However, this result can be extrapolated to 
the thermodynamic limit by a finite-size-scaling analysis [85] [ 86 ]. As shown in Fig. [5](c), 
we find A = 0 for U/t — 0 and A 7 ^ 0 for U/t > 0 in the thermodynamical limit. 
Therefore, here we call the phase, in which n^(k) has sharp edges at k = ±kp and A = 0 
(thermodynamical limit), the metallic phase, whereas the phase, in which n^(k) become 
smoother, all sharp edges disappear, and A ^ 0 (thermodynamical limit), is referred 
as the Mott insulator Since in this Mott insulator the spin orientations of the 
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Figure 7. The spin-correlation functions s(r) for the different SOC strengths (a) 
X/t = 0, (b) X/t = 0.1, (c) X/t = 1, and (d) X/t = 1.73 at the half filling (n = 1). In 
all subfigures, the lattice length is given by L = 100. 


nearest-neighbor sites are antiparallcl, the phase is finally called the antiferromagnetic 
Mott insulator (see also the following discussions). 

In the presence of SOC (see, for example, X/t = 1), the momentum distribution 
function n^{k) is unconventional and A = 0 (thermodynamical limit), when Uft = 0; see 
the black solid curve in Figs. [5](d)-|5](f). It implies that the system is located at the SOC- 
induced metallic phase. When increasing the on-site repulsive interaction strength U/t , 
these momentum distribution functions rq(fc) also become smoother, all sharp edges also 
disappear, and A ^ 0 (thermodynamical limit), i.e., the system enters into the Mott 
insulator. However, as will be shown in the next subsection, these Mott insulators, 
without SOC or with SOC, are quite different. Without SOC, the Mott insulator is 
antiferromagnetic, whereas it becomes spin-rotating antiferromagnetic if 0 < X/t < 1, 
and spin-rotating ferromagnetic if X/t > 1. 

In Figs. E](a)-l6](c), we plot the momentum distribution functions n^(k) and the gap 
A for the different filling factors n, when X/t = 1 and U/t = 10.0. If n ^ 1, all sharp 
edges (Fermi points) in the momentum distribution functions n^{k) still exist and the 
zero gap, A = 0 (thermodynamical limit), remains. These indicate that the system is 
always located at the SOC-induced metallic phase for any on-site repulsive interaction 
strength U/t. However, at the half filling {n — 1), the momentum distribution function 
rq-(fc) becomes smoother and A ^ 0 (thermodynamical limit) for any on-site repulsive 
interaction strength U/t. These mean that these Mott insulators only occur at the half 
filling (n = 1), which is similar to the result without SOC [i5j. 











Spin-orbit coupled repulsive Fermi atoms in a one-dimensional optical lattice 


13 


(a) Xlt= 0.0 

(b) }Jt=0A 

f i V i' o ^ ^ 't r i f i' o \ v -v '/ f i ^ ^ ^ ^ 

(c) XIt= 1.0 

I*— t —»l«— 1 —«l«— I —*1*— I -»l«- 1 -*l*. I -*l«- ^ -*j*- T -»l«- 1 ~*j*~ ^ -p ^ 

(d) A,/t=1 .73 

'\/"\ ^\r' Yr' Yr' ; ~\r' ; ~\Y r f ~ 


Figure 8. Vector plots of DMRG results Si = (sf, sf) for the different SOC strengths 
X/t at the half filling (n = 1). In all subfigures, the on-site repulsive interaction 
strength and the lattice length are given respectively as U/t = 10 and L = 100, which 
are the same values as the red curve of Fig. [7] 


4-2. Spin-correlation function and spin-structure factor 

In Figs. [?la)-[7](d), we plot the spin-correlation functions s(r) for the different SOC 
strengths X/t. These figures show clearly that in the presence of the on-site repulsive 
interaction (U/t ^ 0), quasi-long-range spin correlation appears, i.e., s(r > 1) ^ 0. 
This is contrast to the result without the on-site repulsive interaction, in which only 
the short-range spin correlation emerges; see Figs. |4](a) and S(b). When increasing 
the SOC strength X/t, the spin-correlation functions between the nearest-neighbor 
sites vary from the negative to the positive. In order to see this physics more 
clearly, we plot, in Figs. EXa)-jHXd) , spin distributions of each sites, i.e., si = ( sf,sf ), 
where sf = (G| cLc^ + cj,cq-1 G) and sf = (G| cLq-j- — C4C4 |G) with the ground-state 
wavefunction |G), for the different SOC strengths X/t. We define an intersection angle 
6 between the different spin orientations of the nearest-neighbor sites. 

I 11 the absence of SOC (X/t = 0), the spin orientations of the nearest-neighbor 
sites are antiparallel and 9 = 7 r; see Figs. | 8 ](a). Moreover, the spin-correlation function 
s(r) > 0 if r is odd, while s(r) < 0 if r is even; see Fig. [7|a). These mean that the 
corresponding spin-spin interactions of the nearest-neighbor sites are antiferromagnetic. 
In addition, the spin-correlation function s(r ) decays as a power law and changes the 
sign with a period t = 2 ng ; see also Fig. E(a). This phase is usually called the 
antiferromagnetic Mott insulator. In the presence of SOC (see, for example, X/t = 0.1), 
the spin orientations of the nearest-neighbor sites are not antiparallel (i.e., the spins 
are rotating) and 7 t /2 < 6 < 7 r; see Figs. [ 8 ](b). Since in this case the traditional spin- 
independent hopping still plays a dominate role, the quasi-long-range antiferromagnetic 
spin correlation remains, but with a period 2 < T < 4; see Fig. [71(b). Thus, we 
call the corresponding phase the spin-rotating antiferromagnetic Mott insulator. When 
X/t = 1, the spin orientations of the nearest-neighbor sites are vertical and 6 = 7 t/ 2 ; 
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Quasimomentum k 


Figure 9. (a) The spin-structure factors S(k) for the different SOC strengths X/t 
as functions of the quasimomentum k. (b) The momentum fc pea k as a function of 
the SOC strength X/t. (c) The spin-structure factors S(k) for the different on-site 
repulsive interaction strengths XJ/t as functions of the quasimomentum k, when the 
SOC strength is given by X/t = 1. In (a)-(c), the lattice length is given by L = 100. (d) 
The momentum fc pea k as a function of 1 /L. In (a), (b), and (d), the on-site repulsive 
interaction strength is given by U/t = 10. In all subfigures, the filling factor is given 
by n = 1. 


see Figs. m c). Moreover, s(l) = s(3) = • • • = 0 with a period T — 4, which means that 
no spin correlation between the nearest-neighbor sites can be found; see Fig. Etc). At 
the strong SOC (X/t > 1), the SOC-induced hopping plays a dominate role. In this 
case, the spin orientations of the nearest-neighbor sites tend to parallel and 9 < 7t/2; 
see Figs. E^d). Moreover, s(l) > 0, which indicates that the spin-spin interactions of the 
nearest-neighbors sites become ferromagnetic. In addition, the quasi-long-range spin- 
correlation function s(r ) also decays as a power law, but changes the sign with a period 
T > 4; see Fig. Eld). We call the corresponding phase the spin-rotating ferromagnetic 
Mott insulator. From above discussions, we argue that SOC can drive the system from 
the spin-rotating antiferromagnetic Mott insulator to the spin-rotating ferromagnetic 
Mott insulator. 

Figure E^a) shows the experimentally-measurable spin-structure factors S(k) for 
the different SOC strengths X/t at the half filling (n = 1). In the absence of SOC 
(X/t = 0), the system has the antiferromagnetic order [78l[87|, and the spin-structure 
factor S(k) has a peak at the momentum /c pea k = tt, as expected. When increasing the 
SOC strength X/t, the peak still exists and varies as 

t \JF -T A 2 

kpezk = 2fc max - 7T = 4 arctan(---) - ir. (15) 
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Figure 10. Schematic phase diagram as a function of the SOC strength X/t and the 
on-site repulsive interaction strength U/t at the half filling (n = 1). The black dot 
denotes the metallic phase (MF), the blue line denotes the antiferromagnetic Mott 
insulator (AFMI), and the red line denotes the SOC-induced metallic phase (SMF). In 
the other abbreviations, SAFMI and SFMI denote the spin-rotating antiferromagnetic 
Mott insulator and the spin-rotating ferromagnetic Mott insulator, respectively. 


It is easy to find from Eq. (fT5]l that when X/t = 1, k pea k = 7 t/ 2; see Fig. [9(b). When 
0 < X/t < 1, the system has the spin-rotating antiferromagnetic Mott insulator with 
7r/2 < fcpeak < 7T- When 1 < X/t < 2, the system has the spin-rotating ferromagnetic 
Mott insulator with k pea k < 7 t/ 2. Since the antiferromagnetic order has been detected 
experimentally by measuring k pea k via the time-of-flight imaging [75], these spin-rotating 
ferromagnetic and antiferromagnetic orders can also be detected by the same method. 
For a fixed SOC strength X/t = 1, when increasing the on-site repulsive interaction 
strength U/t, the momentum k pe a k remains unchanged, while the magnitudes of peaks 
increase; see Fig. [9(c). In Fig. [9(d), we present a finite-size-scaling analysis of the 
momentum fc pea k- This figure shows that the momentum k pea k remains unchanged when 
increasing the lattice length L. 

4-3. Phase diagram 

Based on the predicted properties of the momentum distribution, the gap, and the spin 
correlation, we find that the homogeneous Hamiltonian (|3]) has five phases, including the 
metallic phase, the SOC-induced metallic phase, the antiferromagnetic Mott insulator, 
the spin-rotating antiferromagnetic Mott insulator, and the spin-rotating ferromagnetic 
Mott insulator, at the half filling (n = 1). In the metallic phase, the system has 
conventional momentum distribution with the gap A = 0 (thermodynamical limit), 
while it becomes filling-dependent unconventional momentum distribution but with 
the same zero gap (thermodynamical limit) in the SOC-induced metallic phase. In 
the antiferromagnetic, spin-rotating antiferromagnetic, and spin-rotating ferromagnetic 
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Mott insulators, the system has a nonzero gap (thermodynamical limit) and the quasi- 
long-range spin correlation decaying as a power law. However, in the antiferromagnetic 
Mott insulator, 9 = tt, T = 2, and k pea k = tt. In the spin-rotating antiferromagnetic 
Mott insulator, tt/2 <9 <tt, 2 <T <4, and n/2 < k peak < tt. Whereas in the 
spin-rotating ferromagnetic Mott insulator, 9 < tt/2, T > 4, and k pea k < tt/2. 

In Fig. [lOl we give a schematic phase diagram as a function of the SOC strength 
X/t and the on-site repulsive interaction strength U/t. In this figure, the Mott insulator 
and the metallic phase are separated by the gap A (thermodynamical limit) [45]. The 
metallic phase and the SOC-induced metallic phase are separated by the momentum 
distribution function rif(k). The antiferromagnetic Mott insulator and the spin-rotating 
antiferromagnetic Mott insulator are separated by 9 = tt, T = 2, and fc pea k = tt- The 
spin-rotating antiferromagnetic Mott insulator and the spin-rotating ferromagnetic Mott 
insulator are separated by 9 = tt/2 , T = 4, and k pe a k = tt/2. 

We also find that the transitions between the metallic phase and the 
antiferromagnetic Mott insulator, between the SOC-induced metallic phase and the spin- 
rotating antiferromagnetic Mott insulator, and between the SOC-induced metallic phase 
and the spin-rotating ferromagnetic Mott insulator are of second order, because the 
second-order derivative of the ground-state energy is discontinuous at the critical points. 
For the other transitions, the ground-state energy and its derivative vary smoothly. It 
means that no phase transition can be found, although for the different phases, the 
momentum distribution and the spin correlation are different. 

5. Discussions and conclusions 

Before ending up this paper, we make one remark. In the real experiment, the Zeeman 
held usually exists. When the Zeeman held is considered, an extra Hamiltonian 



(16) 


should be added in the Hamiltonian (]3]). This Zeeman held can lead to spin hipping 
in the same site. Moreover, it has a strong competition with SOC, which makes Fermi 
atoms hop between the nearest-neighbor sites with spin hipping, and the conventional 
spin-independent hopping between the nearest-neighbor sites. As a consequence, rich 
magnetic properties can emerge. The deep understanding of relevant behavior is very 
complicated (we need introduce more physical quantities), but interesting. We leave 
this important problem for further investigation. 

In summary, we have investigated spin-orbit coupled repulsive Fermi atoms in a 
ID optical lattice by the DMRG method. We have found that SOC can generate the 
filling-dependent unconventional momentum distribution, whose corresponding phase is 
called the SOC-induced metallic phase, and drive the system from the spin-rotating 
antiferromagnetic Mott insulator to the spin-rotating ferromagnetic Mott insulator. 
We have predicted a second-order quantum phase transition between the spin-rotating 
ferromagnetic Mott insulator and the SOC-induced metallic phase at the strong SOC. 
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Finally, we have also shown that the momentum, at which peak of the spin-structure 
factor appears, can be affected dramatically by SOC. The analytical expression of this 
momentum with respect to the SOC strength has also been derived. Attributed to the 
recent experimental realization of the spin-orbit coupled Bose-Einstein condensates in 
the ID optical lattice [69], we expect that our predictions could be observed in the near 
future. In particular, the spin-rotating ferromagnetic and antiferromagnetic correlations 
can be detected by measuring the SOC-dependent spin-structure factor via the time-of- 
flight imaging m- 
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